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I.  INTRODUCTION 


Previous  investigators  have  studied  the  interactions  of  thin  cylin¬ 
ders  in  free  space  [1],  and  thin  cylinders  over  perfectly  conducting 
ground  [2]  with  an  electromagnetic  step  plane  wave.  This  study  will 
present  in  an  approximate  manner  the  interaction  of  a  thin  cylinder  over 
a  finitely  conducting  ground  with  an  electromagnetic  step  plane  wave. 
Analysis  will  be  conducted  using  the  Singularity  Expansion  Method  (SEM). 

Current  induced  by  the  incident  field  on  the  cylinder  surface  is 
approximated  by  an  axially  directed  filamentary  current  on  the  cylinder 
axis.  Boundary  conditions  are  applied  only  to  the  axial  component  of  the 
incident  electromagnetic  field  on  the  wire  axis.  The  above  "thin-wire" 
assumptions  are  valid  provided  the  length  of  the  cylinder  is  much  greater 
than  its  radius,  and  the  cylinder  is  many  radii  away  from  the  ground 
plane  [3]. 

The  scattered  field  reflected  from  the  ground  plane  is  scaled  by  the 
complex  Fresnel  reflection  coefficient  for  the  appropriate  angle  of  inci¬ 
dence  and  polarization  involved.  Since  the  Fresnel  reflection  coefficient 
is  strictly  valid  only  for  plane  wave  incidence,  the  scattered  field  must 
approximate  a  plane  wave  at  the  free  space-lossy  ground  interface.  In  a 
study  by  Sarkar  and  Strait  [4]  it  was  shown  that  the  above  method,  termed 
"reflection  method"  gave  results  in  the  real  frequency  domain  within  10« 
of  the  exact  Sommerfeld  formulation  for  a  horizontal  electric  dipole  as 
long  as  the  dipole  was  at  least  (0.25x//e^)  from  the  ground  plane.  Speed 

of  computation  is  an  inherent  advantage  when  using  the  reflection  method 
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II.  THEORY 


Integro-Differential  Equation 

The  system  defined  by  Figure  1  is  composed  of  a  thin  cylinder  over 
an  imperfectly  conducting  ground  plane  and  the  incident  electromagnetic 
radiation.  As  shown,  the  cylinder  is  of  length  £,  radius  a,  and  height  n 
above  the  ground  plane.  The  imperfectly  conducting  nonmagnetic,  u  =  uQ, 
ground  plane  is  characterized  by  its  conductivity  a,  and  permittivity, 
e  =  epeQ.  The  incident  plane  wave  is  vertically  polarized  and  propagates 
at  an  angle  9  with  respect  to  the  normal. 

Currents  will  be  induced  on  the  cylinder  by  the  incident  field  and 
by  reflection  of  the  incident  field  from  the  ground  plane.  These  induce! 
^rrents  will  reradiate  in  the  presence  of  the  imperfectly  conducting 
ound  producing  a  primary  and  secondary  scattered  field.  Primary  scat¬ 
tered  radiation  is  that  part  of  the  current  induced  field  which  propa¬ 
gates  directly  from  the  scatterer  surface  to  observation  point.  Secondary 
scattered  radiation  reflects  from  the  ground  plane  before  reaching  the  ob¬ 
servation  point,  and  therefore  must  be  scaled  by  the  complex  Fresnel  re¬ 
flection  coefficient  for  the  appropriate  angle  of  incidence  and  polariza¬ 
tion  involved.  The  total  scattered  field  consist  of  the  sum  of  the  pri¬ 
mary  and  scaled  secondary  parts. 

The  primary  scattered  electric  field  is  related  to  the  currents  that 
produced  it  through  a  magnetic  vector  potential  and  appropriate  differen¬ 
tial  operators.  This  relation  may  be  written  as 


4 


(2.1) 


(2.2) 


y  =  s/c,  c  =  speed  of  light  in  free  space. 

Use  of  the  complex  frequency  variable  s  in  the  above  equations  implies 
that  the  Maxwell  equations  have  been  Laplace  transformed. 

Secondary  scattered  electric  radiation  may  also  be  written  in  terms 
of  a  magnetic  vector  potential  as 

Is  ■=  -vQsP  +  v (v  •  As)/s eQ  (2.3) 

~  secondary  scattered  electric  field 


I 


where 


cv  K. 

-  ■  U/ 


-yIF-F'I5 


K(r*  ,s)  - - r  ds1 

|r-F'  |s 


(2.4) 


=  secondary  magnetic  vector  potential 


_  _  s 

| r- r ' |  =  distance  from  source  point  to 

ground  plane  to  observation  point. 


All  other  terms  in  (2.3)  and  (2.4)  have  been  previously  defined. 

With  surface  currents  on  the  cylinder  approximated  by  a  z-directed 
filamentary  current  on  the  cylinder  axis,  the  primary  and  secondary  mag¬ 
netic  vector  potentials  of  (2.2)  and  (2.4)  reduce  to 


jc 

i-bf 

o 

l 

*-bf 


I (z '  ,s) 


I(z' ,s) 


-y|F-7'|P 


I  r-r ' |p 


■yI r-r' |s 


I  r-r1 1 s 


(2.5) 


(2.6) 


Locate  the  observation  point  on  the  cylinde"  surface,  then 


k-r'|P  -  r,  «  [(z-z')2  +  a2f2 


k-r'|S  =  r9  =  [(z-z')2  +  (2h+a)2f2 


(2.7) 


(2.8) 


An  expression  for  the  primary  scattered  electric  field  in  terms  of  the 
unknown  induced  currents  may  be  formed  by  substituting  (2.5)  into  (2.1), 


the  result  is 


Ep 


bf 


I  (z ' ,s) 


dz'  .  (2.9) 
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A  similar  expression  for  secondary  scattered  electric  radiation  may  be 

formed  by  substituting  (2.2)  into  (2.3)  with  the  result 

£  r 

E!  <z-s>  ■  (-V*ir-4>  .  (2.io) 

0  a  Z  •/  Z 

0 

As  stated,  the  total  scattered  field  consists  of  the  sum  of  the  primary 
and  scaled  secondary  parts,  that  is 


rto t  /  \ 
Ez  (z,s) 


EP(z,s)  -  rJe*(z,s) 


£  -yr 

=  4)fc/I(2''s)1rrdz'  <2-"> 

0  oZ  i 

0 

£  -yr 

-  h>  b  f  I(z'-S>  TT-dz'  • 

0  3Z  ./  2 


The  complex  Fresnel  reflection  coefficient  for  vertical  polarization, 
i.e.  (Rp  is  used  in  (2.11)  to  scale  the  z-directed  secondary  electric 
field,  and  is  defined  as 


(er+X)  sinij;  -  [(eR+X)  -  cos^ifi]^ 
(eR+X)  s i nip  +  [ ( e  +X)  -  cos^V] 2 


(2.12) 


where  X  =  120Tra/y,  y  =  s/c. 

This  expression  utilizes  Jordon's  [5l  definition  of  the  plane  wave 
reflection  coefficient  for  vertical  polarization.  As  pointed  out  by  a 
previous  investigator  [6],  the  minus  sign  n  (2.11)  comes  from  Jordon's 
assumed  positive  directions  of  electric  fields  for  the  incident  and  reflec¬ 
ted  waves  [5].  Note  that  in  (2.12)  is  the  angle  formed  by  the 


i 
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secondary  scattered  incident  ray  and  the  ground  plane.  It  is  a  function 
of  source  and  field  point  position  as  depicted  by  Figure  2-1. 

In  Figurj  2-2  is  shown  the  vertically  polarized  total  incident  and 
transmitted  electromagnetic  excitation.  At  the  free-space  lossy  ground 
interface  the  tangential  components  of  the  total  incident  field  must 
equal  the  tangential  components  of  the  transmitted  field.  This  condition 
is  equivalent  to  requiring  normal  wave  impedances  to  be  continuous  at  the 
interface  [7],  thus 


where 


Z^  =  x-directed  wave  impedance  in  free  space  region 

•  _•  -y (-x  sine .+z  cose . ) 

=  E^s)  sinei  [e  1  1 

-y(x  sine .  +  z  cose . ) 

-re  1  1  ] 


(2.14) 


=  total  z-directed  incident  electric  field 

_•  ^  -y(-x  sine.+z  cose.) 

Hy  •  Ej(s)/ni  [e 

-y(x  sine.  +  z  cose.) 

+  r  e  ’  '1  (2., 5) 


=  total  y-directed  incident  magnetic  field 


=  x-directed  wave  impedance  in  lossy  region 
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Ez  =  E^s)  sinet  e 


-y2(-x  sinet+z  coset) 


(2.16) 


=  z-directed  transmitted  electric  field 

.  -Yo(-x  sine.+z  cose. ) 

-  E^(s)/n2  e  2  t  1 


=  y-directed  transmitted  magnetic  field 


(2.17) 


n-j  =  intrinsic  impedance  of  free  space 

n2  =  intrinsic  impedance  of  conducting  earth 

Y2  =  complex  propagation  number  of  conducting  earth 

and  r  is  the  reflection  coefficient  to  be  determined.  Note  that  the  total 
incident  field  consists  of  the  direct  electromagnetic  excitation  plus  its 
reflection  from  the  ground  plane.  The  principal  of  direction  cosines  has 
been  used  in  the  development  of  (2.14)  thrcugh  (2.17).  Substituting  the 
field  expressions  into  (2.13),  and  simplifying  the  result,  gives 


-Y2hsine . 

r  =  e  1 


n-j  sinei  -  n2  sinet 
~ni  sinei  +  ~n2  "sinFJT 


(2.18) 


Through  the  application  of  Snell's  law  of  refraction  (2.18)  may  be  written 
-Y2hsine  (eD+X)  sine  -  [(eD+X )-cos^e]^ 

Y  S  0  _ _ _ _ _ 

(eR+X)  sine  +  [(eR^)-cos2er2 

-Y2hsine 

R! 


(2.19) 
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Equation  (2.19)  differs  from  (2.12)  only  by  the  exponential  factor  and 
angular  dependence  of  R  .  The  total  z-directed  incident  electric  field 
may  now  be  obtained  by  substituting  (2.19)  into  (2.14)  with  the  result 

.  .  -y(-x  sine+z  cose) 

e’  =  E^(s)  sine  [e  -  rJ 

-y(2h  sine+x  sine+z  cose) 

e  ]  .  (2.20) 

To  insure  uniqueness  the  total  z-directed  scattered  electric  field 
must  cancel  the  total  z-directed  incident  electric  field  on  the  cylinder 
axis,  therefore 

EzJ  =  "  Ez°t(z>s)J  ’  (2,21) 

x=o  x=o 


Evaluating  tha  z-directed  scattered  and  incident  fields  on  the  cylinder 

axis  rather  tnan  on  the  surface  is  an  appropriate  thin  wire  assumption. 

Using  (2.11)  and  (2.20)  in  (2.21)  and  applying  the  above  discussed  thin 

wire  assumption,  (2.21)  becomes 

.  -yz  cose  -y (z  cose+2hsine) 

eJ(s)  sine  [e  -  R®  e  ] 


£  -  yT 

(wos"  st  7J]  bf  I(z’’s)  ~T7~dz' 

Ode 


l 


-Yrc 


r:'(v"5 k^'^f  I(2',s) 


(2.22) 


Define 


-yr. 


Fp(z,z',s)  = 


=  primary  Green's  function 


(2.23) 


11 


-yr2 

Fs(z,z',s)  =  - - =  secondary  Green's  function  (2.24) 

r2 


where 

r,  =  t(z-z')2  +  a  2fz 
r2  =  [(z-z1)2  +  (2h+a  )2f> 

and  since 


1  (-4’V  <  v  -  d-  ^  ■  v2) 


4-n 


(2.25) 


the  integro-differential  equation  to  be  solved  for  the  unknown  induced 
currents  is 


~yz  cose  -y(z  cose+2h  sine) 

(-4ire  s)  E  (s)  sine  [e  -Re 

o  o  v 


] 


t 

dz‘ 


,  2, 
( — 2  "  Y  ' 


£ 

/ 


I (z ' ,  s )  Fs (z,z 1 ,s)  dz1 


2  > 

-rJ(^-o-y2)  /  I(z\s)  RS(z,z',s)  dz'  (2.26) 

92  Jo 

Application  of  the  Method  of  Moments 
The  purpose  of  this  section  is  to  reduce  the  integro-differential 
equation,  (2.26),  to  a  form  suitable  for  numerical  solution.  The  tech¬ 
nique  by  which  this  may  be  accomplished  is  known  as  the  method  of  mo¬ 
ments  [8]. 

Generally,  the  method  of  moments  may  be  used  to  solve  an  inhomoge¬ 
neous  equation 


L(f)  =  g 


(2.27) 
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where  L  is  a  linear  operator,  g  is  known,  and  f  is  to  be  determined.  The 
unknown  f  is  expanded  in  a  series  of  functions  as 


(2.28) 


where  the  f  are  called  basis  functions  and  the  a  are  constants.  Using 
n  n  3 

the  linearity  of  L  and  (2.28),  (2.27)  may  be  written 


E«„  W„>  ■  9 


(2.29) 


Next  a  set  of  testing  functions  is  defined,  w^ ,  w^,  w^  ...,  and  the  inner 
product  of  (2.29)  with  each  wm  formed,  yielding 

£w  L(fn)>  =  <wm’9>  •  (2-30) 

n 


n  =  1 ,  2,  3,  ... 
m  =  1 ,  2,  3,  . . . 

Equation  (2.30)  is  a  matrix  equation  which  may  be  solved  using  standard 
numerical  techniques. 

For  the  problem  at  hand,  a  geometric  interpretation  of  the  method 
of  moments  may  be  forwarded.  Let  the  thin  wire  be  broken  into  segments, 
with  the  induced  current  on  each  segment  assumed  to  be  constant.  See  Fig¬ 
ure  2-3.  This  is  equivalent  to  expanding  "he  current  in  a  set  of  pulse 
functions  as 

1  (z 1  ,s )  «n(s>  !„(*■)  (2.31) 

n 


14 


where 


a  (s)  =  unknown  coefficient  of  constant  current 

n 


Uz')  = 


in  the  n  subsection, 
(l  for  zn  <  z  <  zn+^ 


0  elsewhere  n  =  2,  3,  ...  N-l 


(2.32) 


At  wire  ends  the  boundary  condition  of  zero  current  may  be  satisfied  by 
defining 


I(o)  =  I(£)  =  0 


(2.33) 


Actually,  this  condition  is  automatically  met  by  allowing  the  two  end 
subsections  to  extend  past  the  end  of  the  scatterer.  Using  a  pulse 
function  expansion  of  the  unknown  induced  current  allows  the  integration 
in  (2.26)  to  be  approximated  by  a  sum  of  integrations  over  N  segments. 


Thus, 


.  -yz  cose  -y(z  cose+2h  sine) 

(-4teqs)  E^(s)  sine  [e  "  Rv  e  ■ 


E“n(s)  f 
II  A 


(— o  "  Y2)  Fp(z,z',s)  dz ' 
dz 


(2.34) 


n  J  n 


(■^-j  -  y^)  Fs(z,z',s)  dz1 

3Z 


where 


A  =  £/(N-l)  =  length  of  a  zone 


(2.35) 


i 
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N  =  number  of  subsections  or  zones 


zn  =  (n  -  3/2 ) A.  n  =  1,  2,  ...  N+l 


=  subsection  ends. 


A  set  of  testing  functions  may  now  be  defined  as 

1  z=z 

6(2-zm)  - 

0  ztz„ 


(2.36) 


where  zm  =  (m-l)&  m  *  1,  2,  ...  N 


(2.37) 


Now  by  forming  the  inner  product  of  (2.34)  with  each  of  the  delta  func¬ 
tions  of  (2.36),  and  approximating  the  differentiation  in  (2.34)  by 
finite  differences,  that  is 


-  — ]—p  [f ( z+az )  -  2F(z)  +  F(z-Az)] 
dz2  (Az  j 

the  integro-differential  equation  becomes 


(2.38) 


,  -yz  cose  -y(z  cose+2h  sine) 

(-4,e0s)  e’(s)  sine  [e  m  -  R®  e 


E  “n(S>  7  «  f  FP( 

n  A  <n 


zm+1,z',s)  -  (y2A2+2)Fp(zm,z’ ,s) 


+  Fp(zm_rz',s)  dz1]  -  Rj  [  f 

J  n 


F  <WZ'-S> 


-(y242+2)  Fs(zm,2',s)  *  Fs(zml,z' ,s)  dz']) 

n  =  2,  3,  ...  N-l 
m  =  2,  3,  ...  N-l . 


(2.39) 
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Equation  (2.39)  is  a  set  of  linear  algebraic  equations,  and  may  be 
placed  in  the  form  of  a  matrix  equation, 

V(s)  =  f(s)  T(s),  (2.40) 


where  a  single  bar  represents  a  column  matrix  or  vector  and  double  bars 
indicate  a  square  matrix.  Define  the  matrices  as 


where 


V(s)  =  the  source  vector  =  [v  ], 


vm  ~  the  matrix  elements  of  V(s) 

,  -yzm  cose  -y(z  cose+2h  sine) 

=  (-4,cos)  E^(s)  sine  [e  m  -  rJ  e  m 


] 


(2.41) 


m  =  2,  3,  ...  N-l 


I (s )  =  the  response  vector  =  [i'n] 


where 


ip  =  the  matrix  elements  of  I(s) 

=  an>  unknown  coefficients  of  constant 

current  in  the  nth  zone  (2.42) 

n  =  2,  3,  ...  N-l; 

Z(s)  =  the  impedance  matrix  =  [zn,nl> 
where  zm  =  the  matrix  elements  of  7(s) 
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zmn  A° 
A 


zn+1 

?"/ 

•'ll 


fP(V1’z,*s)  -  (y2A2+2)  Fp(zm,z\s) 


.n+l 


+FP(2m.1.zl.s)dz,]-Rl[ 


L 


F  (Wz'-s) 


(T2a2*J)  Fs(zm,z',s)  +  Fs(zb1_],z',s)  <Jz']} 
n  =  2,  3,  ...  N-l 


m  =  2,  3,  ...  N-l 
The  unwieldy  appearance  of  (2.43)  may  be  improved  by  defining  two 
functions  n+l 


(2.43) 


Hp  n(z  , s)  =  f 
m,n  m  '  I 

Jzn 

Hm,n<Vs)  =  [ 
J_  n 


Fp(zm,z\s)  dz' 


F  (z.z'.s)  dz1 


-y[(z-z')2  +  a  2f 


with 


m 


FH(z  , z' ,s)  =  —  _ 

m  r  /  ,  ?  ?  J- 

[(zm-z')2  +  a2]2 


F  (VZ'’S) 


-y[(z  -z')2  +  (2h+a )2f2 


[(z  -z')z  +  (2h+a )^] 


(2.44) 


(2.45) 


(2.46) 


(2.47) 


mn 


Now  (2.43)  may  be  redefined  as 

?  (ln<Ws>  -  (Y2A2+2)H^>n<zm.s) 

+  OVl’s>  *  K  C>£.„<WS>  -  (T2A2*2)HF_n(zm,s) 


HS  (z  ltS)]}  n  =  2,  3,  ...  N-l 
m»n  m-l’s;j>  m  =  2,  3,  ...  N-l 


(2.48) 
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Alternate  Formulation  for  the  Secondary  Scattered  Electric  Field 

As  pointed  out  by  Shumpert  [9],  an  alternate  formulation  for  the 
secondary  scattered  electric  field  is  derivable  by  using  the  exact 
expression  for  the  electric  field  produced  by  a  constant  current  element. 
Using  (2.11)  in  (2.21)  gives 

-E l  =  EP(z,s)  -  rJe^z.s)  (2.49) 

where  eP(z,s)  and  R^E^(z,s)  are  the  primary  and  scaled  secondary  scattered 
electric  fields  on  the  cylinder  surface  respectively,  due  to  the  unknown 
induced  current.  Since  the  wire  scatterer  has  been  segmented  as  shown 
in  Figure  2-3,  and  the  current  induced  on  each  segment  assumed  constant, 
it  follows  that  the  secondary  scattered  electric  field  will  be  merely 
the  sum  of  the  fields  caused  by  each  individual  constant  current  element. 
The  exact  electric  field  due  to  a  constant  current  element  has  been  given 
by  Harrington  [7].  Applying  the  principles  discussed  above,  and  using 
Harrington's  expression  for  a  constant  current  element,  the  secondary 
scattered  electric  field  may  be  written 


Ez(z,s) 


V 

2tt 


N-l 

L 

n=2 


-yr 


mn 


t(7- 

mn 


2 

COS  \f>, 


Yr, 


mn 


mn 


r(~t— 


mn 


1 

rL 

mn 


yr 


2 

sin  ] 
mn 


mn 


(2.50) 
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The  obvious  utility  of  the  alternate  expression  for  zmn  is  that  an  inte¬ 
gration  has  been  replaced  by  an  algebraic  expression,  and  thus  numerical 
evaluation  tine  reduced. 

Application  of  the  Singularity  Expansion  Method 
The  Singularity  Expansion  Method  (SEM)  formalized  by  Baum  [10-11], 
and  applied  by  many  others,  allows  one  to  treat  a  distributed  system, 
such  as  the  one  at  hand,  in  a  manner  similar  to  that  used  in  classical 
circuit  theory.  In  circuit  theory  the  time  domain  response  of  a  linear 
circuit  excited  by  an  arbitrary  waveform  may  be  determined  by  knowledge 
of  the  location  of  any  singularities  of  the  response  function  as  well  as 
the  corresponding  residues  [1].  In  the  case  of  a  distributed  system 
there  are  an  infinite  number  of  singularities,  and  associated  with  each 
is  a  natural  modal  current  distribution.  For  any  arbitrary  excitation, 
i.e.,  incident  electromagnetic  radiation,  one  need  only  determine  how 
much  of  each  natural  modal  current  has  been  excited  [12].  This  is  deter¬ 
mined  by  the  coupling  coefficient  associated  with  the  given  singularity. 
The  solution  of  (2.40)  is 

T(s)  =  Z_1(s)  V(s)  =  Y(s)  V(s)  (2.53) 

Y(s)  =  the  inverse  of  the  system 
impedance  matrix 


y _  =  the  matrix  elements  of  Y(s). 

•'mn 


where 
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Elements  of  the  inverse  matrix  may  be  written 

hn  A 

nn 
a(s) 


<-’)ra+n  v,<s> 


'mn 


(2.54) 


th 


where  Anm ( s )  is  a  minor  determinant  of  Z(s)  formed  by  deleting  the  n 
row  and  the  m^  column,  and  A(s)  is  the  determinant  of  Z(s).  Examination 
of  (2.54)  indicates  that  the  poles  of  the  response  function  I(s),  are 
the  zeros  of  a(s).  The  poles  have  been  termed  the  system  natural  reso¬ 
nant  frequencies.  Now  writing  the  unknown  current  in  a  partial  fraction 
expansion,  one  obtains 

fr 

T(s)  =  Y(s)  V(s)  J?Ti77  V(s)  .  (2.55) 


t  h  —  *y< 

In  (2.55)  s^  is  the  i  natural  resonant  frequency  and  Y.  the  correspond 
ing  residue,  defined  as 

ttT 


Yi  =  lim  [(s-s.)  Y(s )] 


(2.56) 


s-Si 


Using  the  circuit  theory  analogy,  certain  information  about  ^be  system 
natural  resonant  frequencies  may  be  inferred.  First,  the  r  ances 
must  occur  in  the  left-hand  portion  of  the  complex  plane  to  isure  a 
decaying  response.  The  poles  must  occur  in  conjugate  pairs  to  produce  a 
real  time  domain  current,  and  since  the  scatterer  has  a  finite  quality 
factor,  no  poles  may  reside  on  the  jw  axis.  Furthermore,  it  is  assumed, 
but  without  proof,  that  the  poles  are  all  simple  [1].  This  has  been 
found  to  be  tie  case  in  many  exactly  solvable  geometries. 

The  system  residue  matrix  at  the  pole  s=sa,  Y1",  has  been  shown  to 


be  a  dyadic  [13],  [14];  that  is,  it  is  proportional  to  the  outer  product 
of  two  n  dimensional  vectors 


where 


and 


r  =  b  v  cc 

a  ana 


v  =  natural  mode  vector 

a 

C*  =  transpose  of  coupling  vector 


(2.57) 


=  proportionality  constant. 

The  natural  mode  vector  is  a  solution  to  the  equation 

=  0  (2.58) 

and  the  coupling  vector  satisfies  the  equation 

^sa)  ^  =  0  (2.59) 

where 

Zt(sa)  =  transpose  of  Z(sa)  . 

In  this  problem  the  system  impedance  matrix  is  symmetric,  that  is 

=  Z(sa)  (2.60) 

so  the  natural  mode  vector  and  the  coupling  vector  are  identical, 

(2.61 ) 


and  therefore,  (2.57)  may  be  written 


I 
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where 


v  =  transpose  of  v 
a  a 


Let  the  natural  mode  vector  be  normalized  such  that  its  maximum  element 
is  real  and  equal  to  unity. 


v  =  v  normalized 
“o  a 


(2.63) 


Now  the  residue  matrix  becomes 


a  a  a  a 

0  0 


(2.64) 


Note  that  the  proportionality  constant  0^  is  not  the  same  as  in  (2.62), 
but  no  confusion  should  result  since  it  has  not  yet  been  defined. 

Several  methods  are  available  for  calculating  the  proportional i ty 
constant  Ba.  The  method  presented  here  is  simple  and  requires  a  minimum 
of  computation  time.  From  (2.64)  it  is  evident  that 


(yr)ii  =  3  {(v  )•>  {(v  ).) 

a  11  a  a  i  o  a  l  0 

,2 


•  «Vij; 


(2.65) 


where 


( y^ ) i i  =  an  element  from  the  ith  row  and  ith 


column  of  the  system  residue  matrix 
evaluated  at  the  singularity  sa 
{(vo)i)  =  an  element  from  the  itfl  row  of  the 
natural  mode  vector  corresponding  to 
the  singularity  s 


*pJ"‘  '' 


.A 
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and  therefore, 


. 


(2.66) 


t  h 

The  ij  element  of  the  system  residue  matrix  at  the  pole  s^  may  be 
found  using  (2.54)  and  (2.56) 


(yr ) . .  =  1 im  (s-s  ) (y  ) . . 

VJV  ij  s^s  aMJVij 

a 


=  1  im 
s->s 


<»„> 

X{s) 


(2.67) 


so  that  finally. 


3  =  lim 

“  S-S  {(v).}^A(s) 
a  a  1  0 


(2.68) 


Expressing  the  residue  matrix  as  in  (2.64)  allows  one  to  write  the  re¬ 
sponse  vector  of  (2.55)  in  the  form 


I(s)  =  Y''  6  v  vt  ^ 

/  -j  a  a  oi  (S-S  ) 
0  0  a 

a 


(2.69) 


or  equivalently  as 


I(s)  T  3  v  C1 

L-J  a  a_  a 


-  Qt  V (s ) 

0  «o  <IrV 


(2.70) 


since  by  (2. 61 ) 


v  =  C  =  normalized  coupling  vactor 
a  a 
0  0 


In  (2.70)  the  scalar  product  of  the  normalized  coupling  vector  C*  ,  and 

ao 

the  incident  field  vector  V(s)  occurs,  mul  i pi i ed  by  the  proportionality 
constant  6  .  This  quantity  has  been  defined  by  Baum  [14],  as  the  cou¬ 
pling  coefficient  c^  at  the  singularity  s=sa,  so 


c 

a 


b  r  [v  (s  )] 

a  a  0  a 
0 


(2.71) 


As  stated,  the  coupling  coefficient  determines  how  much  of  each  natural 

modal  current  distribution,  (i.e.  natural  node  vector)  will  be  used  in 

calculating  tne  response  vector  I(s).  It  is  a  function  of  the  incident 

angle  0,  ground  parameters,  (conductivity  a,  and  permittivity,  e^)  and 

the  particular  singularity  at  which  it  is  calculated.  The  response 

vector  may  now  be  written  in  terms  of  the  coupling  coefficient  as 

c  v 

!(s)  t^tt  •  (2-72) 

a 

a 

Consideration  of  (2.72)  will  show  that  the  complex  natural  frequencies, 
and  natural  made  vectors  are  not  a  function  of  the  angle  of  incidence  of 
the  electromagnetic  excitation,  only  the  coupling  coefficient.  Therefor 
once  these  quantities  have  been  determined  for  a  particular  geometry  the 
response  function  for  any  incident  field  is  easily  found. 

Time  domain  currents  may  be  found  by  using  the  Laplace  inversion 


formula,  thus 


(2.73) 


where  is  the  Bromwich  contour. 
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Approximations  and  Limitations  Imposed 

The  validity  of  any  analysis  is  a  function  of  the  approximations 
used  in  its  construction.  This  section  will  discuss  the  limitations 
imposed  by  the  approximations  used  in  this  analysis. 

The  assumptions  and  approximations  are  as  follows: 

a)  Current  is  assumed  to  flow  only  in  the  direction  of  the  wire 
axis. 

b)  Boundary  conditions  are  applied  only  to  the  axial  component  of 
the  electromagnetic  field. 

c)  The  surface  current  density  is  approximated  by  a  filament  of 
current  on  the  wire  axis. 

d)  End  caps  on  the  cylinder  are  ignored. 

e)  The  moment  method  is  an  approximate  numerical  solution. 

f)  The  reflection  method  is  an  approximate  technique. 

For  a  finite  length  cylinder  the  axially  directed  incident  field 
excites  both  an  axially  and  circumferentially  directly  current  [15-18]. 
However,  the  axial  component  of  the  current  is  much  more  significant  than 
the  circumferential  component  provided  the  length  of  the  cylinder  is  much 
greater  than  its  radius.  Therefore,  the  first  two  approximations  are 
valid  for  thin  cylinders.  Replacing  the  surface  currents  with  filamentary 
currents  on  the  wire  axis  is  valid  provided  the  circumferential  variation 
of  the  surface  currents  is  uniform  [3].  This  will  be  the  case  if  the 
cylinder  is  thin  and  located  many  radii  away  from  the  ground  plane.  The 
scattered  field  contributed  by  currents  induced  on  the  ends  of  the  cylin¬ 
der  will  be  regligible  provided  the  cross-sectional  area  of  the  cylinder 
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is  small,  a<<A.  In  applying  the  moment  method,  the  number  of  zones  into 
which  the  scatterer  is  divided  must  be  increased  as  the  frequency  of 
analysis  increases.  For  an  acceptable  solution,  ten  zones  per  half  wave¬ 
length  may  be  used.  Thus,  one  may  conclude  that  this  analysis  is  not 
applicable  for  high  frequencies. 

The  term  "ground  wave"  applies  to  energy  propagated  over  paths  near 
the  earth's  surface  [5].  It  is  convenient  to  divide  the  ground  wave 
into  a  "space  wave"  and  "surface  wave".  The  space  wave  is  made  up  of 
direct  and  ground-reflected  energy.  The  surface  wave  is  that  energy 
which  is  guided  along  the  earth's  surface,  in  much  the  same  manner  as  an 
electromagnetic  wave  is  guided  by  a  transmission  line. 

Sommerfeld  [19],  was  the  first  to  treat  radiation  from  a  vertical 
dipole  over  a  finitely  conducting  earth.  In  this  original  discussion, 
Sommerfeld  stited  that  it  was  possible  to  divide  the  ground  wave  strength 
into  two  parts,  a  space  wave  and  a  surface  wave.  Norton  [20],  later 
expressed  the  fields  for  an  electric  dipole  above  a  finitely  conducting 
earth  in  a  fo~m  which  clearly  showed  this  separation  into  space  and  sur¬ 
face  waves.  As  pointed  out  by  Jordon  [5],  when  the  dipole  is  located  fa1* 
from  the  earth  the  space  wave  becomes  the  total  ground  wave.  But  as  the 
dipole  nears  the  earth  additional  terms  must  be  taken  into  account  in 
order  to  form  the  total  reflected  field.  These  terms  are  the  ones  which 
account  for  the  surface  wave. 

A  comparative  numerical  study  of  several  methods  for  analyzing  a 
vertical  thin  wire  antenna  over  a  finitely  conducting  ground  plane  was 
done  in  a  dissertation  by  Jerry  McCannon  [21].  In  this  work  he 
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reflection  method  was  found  to  give  answers  within  1  to  2  percent  of  the 
exact  solution  when  the  height  of  the  dipole  was  greater  than  3/8  X. 

As  statea  in  the  introduction,  it  was  found  that  the  reflection 
method  for  the  horizontal  dipole  produced  results  within  10^  of  the  exact 
Sommerfeld  formulation  if 

h  >  (0. 25//TT)x  (2.73) 

K 

where 


h  =  height  of  the  dipole  above  ground 
eR  =  relative  permittivity  of  ground 
X  =  free  space  wavelength. 

For  an  average  relative  permittivity  of  15  the  dipole  must  be  located 
at  a  height  greater  than  h  =  .065X. 

In  short,  the  reflection  method  used  in  this  analysis  gives  good 
results  provided  the  scatterer  is  not  brought  into  close  proximity  of  the 
ground  plane. 


III.  Numerical  Results 


A  computer  code  has  been  written  to  determine  specific  SEM  para¬ 
meters  (i.e.,  the  system  resonances  and  corresponding  natural  modal 
current  distributions).  These  parameters  are  a  function  of  scatterer 
height-to-length  ratio,  length-to-radius  ratio,  ground  conductivity  and 
relative  permittivity. 

Before  presenting  the  results  of  this  work,  a  brief  review  of  two 
other  closely  related  problems  is  in  order. 

Tesche  [1]  treated  the  case  of  a  cylindrical  scatterer  in  free 
space.  The  exterior  natural  resonances  of  this  problem  occur  in  layers 
in  the  complex  plane  and  may  be  described  by  yp  where  "l"  denotes  the 
layer  of  the  Dole  and  "n"  the  pole  within  the  layer.  These  free  space 
resonances  are  repeated  in  Figure  3-1.  Singularities  located  with  an 
"x"  have  natural  modal  current  distributions  which  are  even  functions 
about  the  scatterer  midpoint;  those  located  with  a  dot  have  odd  modal 
distributions.  Note  that  the  imaginary  part,  u>£/c,  of  the  singularities 
in  the  first  layer  occurs  at  approximately  an  integer  multiple  of  n,  or 
at  a  point  where  the  length  of  the  cylinder  is  resonant  (i.e.,  l=n\/ 2). 
First  layer  resonances  are  of  greatest  importance  in  calculating  induced 
currents  since  their  position  in  the  complex  plane  is  nearest  the  imagi¬ 
nary  axis.  In  the  free  space  problem  pole  locations  are  a  function  only 
of  cylinder  length-to-radius  ratio.  As  the  radius  of  the  structure  is 
increased,  the  poles  move  away  from  the  imaginary  axis  indicating  that 
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more  damping  is  introduced.  A  commonly  used  shape  parameter  is  defined 

by 

a  =  21 n  {1/ a)  (3.1) 

All  the  figures  in  this  work  are  for  a  shape  parameter  of  10.6  (£/a=200). 

The  problem  of  a  thin  cylinder  over  a  perfectly  conducting  ground 
plane  has  been  treated  by  Shumpert  [9].  Singularities  in  this  problem 
are  a  function  of  scatterer  height-to-length  ratio  as  well  as  length-to- 
radius  ratio.  When  the  cylinder  is  half  its  length  above  the  perfectly 
conducting  ground  plane,  the  poles  will  be  oriented  as  shown  in  Figure  3-2. 
Comparison  of  this  figure  with  Figure  3-1  indicates  that  the  critical 
resonances,  those  along  the  imaginary  axis,  are  only  slightly  displaced 
from  their  free  space  counterpart.  Figure  3-3  displays  the  movement 
of  singularity  as  the  scatterer  recedes  from  the  ground  plane.  As 
shown,  this  pole  spirals  about  the  free  space  location  until  a  pole  from 
another  layer  takes  its  place.  That  is,  tie  original  pole,  y^,  leaves 
the  spiral  path  and  begins  to  approach  the  origin  while  the  new  pole 
takes  up  the  spiral  trajectory  left  by  the  original  pole.  This  pole 
makes  only  a  partial  revolution  about  the  free  space  location  before  it, 
too,  is  replaced  by  a  new  pole.  It  is  interesting  to  note  that  each  of 
these  singularities  have  similar  modal  current  distributions.  First 
layer  singularities  are  associated  with  the  length  of  the  scatterer,  as 
in  the  free  space  case.  Singularities  in  the  other  layers  have  been 
associated  with  scatterer-ground  plane  interactions. 

Let  us  now  consider  the  problem  of  a  thin  cylinder  over  a  finitely 
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conducting  ground  plane.  Figure  3-4  shows  the  movement  of  singularity 
Y-|-|  in  the  complex  plane  as  several  system  parameters  are  varied.  The 
outer  dashed  spiral  through  the  points  labjled  "A"  is  the  trajectory 
of  pole  y.j  -j  as  the  scatterer  is  brought  near  the  ground  plane;  the  con- 

p 

ductivity  is  held  constant  at  o  =  1.2  x  10  along  this  curve.  A 
conductivity  of  this  value  corresponds  to  a  very  good  earth,  that  is, 
the  ground  plane  may  be  considered  perfectly  conducting.  Therefore, 
this  spiral  is  seen  to  be  identical  with  the  trajectory  of  y^  presented 

in  Figure  3-3.  For  each  value  of  h /£,  the  conductivity  is  varied  from 

2  -4 

a  =  1.2  x  10  ,  point  "A"  on  the  dashed  curve,  to  a  =  1.2  x  10  .point 

"G",  intermediate  values  are  shown  at  points  "B"  through  "F".  The  value 

o  =  1.2  x  10~3  and  o  =  1.2  x  10"^  at  points  "E"  and  "F"  respectively 

correspond  to  typical  values  of  conductivity  for  normal  terrain.  Let 

the  paths  traversed  by  the  pole  for  a  given  value  of  h /£  be  called 

the  inner  spirals.  It  is  seen  that  each  inner  spiral,  corresponding 

to  a  given  value  of  h/£,  converges  to  point  "G"  as  the  conductivity  of 

the  ground  is  reduced,  this  result  is  to  be  expected  since  "G"  is  the 

location  of  y^  for  free  space  conditions.  When  the  pole  is  displaced 

from  position  "G"  along  one  of  the  in^r  spirals  more  energy  is  being 

reflected  from  the  ground  plane.  As  a  becomes  very  large,  points  "A", 

all  the  incident  energy  is  reflected,  and  the  problem  becomes  that  of  a 

cylindrical  scatterer  over  a  perfectly  conducting  ground  plane.  Figure 

3-5  is  also  a  plot  of  y^  as  the  value  of  h/£  and  o  vary,  the  relative 

permittivity  is  held  at  five.  Unlike  Figure  3-4,  point  "G"  is  not  the 

same  in  the  limiting  case  of  small  a,  but  rather  each  inner  spiral 
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converges  to  some  point  along  an  inner  dashed  spiral.  Although  the 
conductivity  becomes  small,  the  ground  plane  is  still  somewhat  reflec¬ 
tive  since  tie  relative  permittivity  is  held  at  five.  Next  consider 
Figures  3-6  and  3-7  where  the  relative  permittivity  is  held  at  fifteen 
and  one  hundred,  respectively.  Again,  as  the  conductivity  is  reduced, 
the  paths  followed  by  the  singularity  for  a  given  value  of  h /l  converges 
to  an  inner  spiral.  The  points  along  these  inner  spirals  are  displaced 
from  the  free  space  pole  position  in  proportion  to  the  relative  permit¬ 
tivity  of  the  ground. 

Some  insight  into  the  behavior  of  the  system  singularities  may  be 
obtained  by  considering  their  origin.  As  stated,  the  singularities  in 
the  complex  plane  are  the  zeros  of  the  determinant  of  the  system  imped¬ 
ance  matrix.  Therefore,  as  the  impedance  elements  vary  so  will  the 
position  of  a  given  pole.  The  impedance  elements  for  this  problem  may 
be  expressed  as 
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where 

z^r  =  matrix  elements  of  the  free  space  problem 

z(J]n  =  matrix  elements  due  to  image  terms  in  the 
perfect  ground  problem 


and 


(eR+X)  sim|j  -  [(eR+X)  -  cos^]"5 
(eR+X)  siniii  +  [(eR+X)  -  cos^]^ 


(3.2) 
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ep  =  relative  permittivity  of  earth 

X  =  120tto/y,  a  =  conductivity  of  earth. 

From  (3.2)  it  is  apparent  that  when  the  reflection  coefficient  is  very 
small  the  impedance  elements  will  be  the  elements  of  the  free  space 
problem,  and  thus  the  singularity  locations  will  be  those  of  Figure  3-1. 
This  condition  will  be  met  when  the  relati/e  permittivity  of  the  earth 
is  unity  and  the  term  "X"  in  (3.3)  is  much  less  than  one.  When  the 
Fresnel  reflection  coefficient  is  near  unity  the  impedance  elements  will 
be  those  of  the  perfect  ground  case,  and  the  singularity  locations  will 
be  those  of  Figure  3-2.  This  occurs  when  either  the  term  "X"  or  the 
relative  permittivity  of  the  earth  is  large.  Note  that  the  size  of  the 
term  "X"  is  proportional  to  the  conductivity  of  the  earth  and  inversely 
proportional  to  the  frequency  of  the  particular  singularity  under 
consideration.  Thus  for  higher  order  poles  the  ground  conductivity 
must  be  larger  to  produce  a  perfectly  conducting  earth  than  for  lower 
order  poles.  In  short,  displacement  of  a  singularity  from  its  free 
space  position  is  a  function  of  the  magnitude  of  discontinuity  in  the 
ground  plane  whether  it  be  produced  by  the  conductivity  or  relative 
permittivity. 

Trajectories  of  first  layer  singularities  y-|4»  and  y^$ 

are  presented  in  Figures  3-8  through  3-15.  As  the  system  parameters  vary, 
these  singularities  behave  similarly  to  the  fundamental  resonance  dis¬ 
cussed  above. 

Figures  3-16  through  3-26  were  constructed  in  order  to  determine  the 
percent  change  in  the  fundamental  damping  constant,  (ReY-j-|)»  relative  to 
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the  free  space  damping  constant,  (R  yQ ) ,  and  the  percent  change  in  the 
fundamental  resonant  frequency,  (Imy^),  relative  to  the  free  space 
resonant  frequency,  (Imy  ).  This  data,  talen  from  the  trajectories  of 
Figures  3-4  through  3-11  is  presented  as  a  function  of  conductivity  with 
the  relative  permittivity  and  height  of  the  scatterer  above  the  ground 
plane  as  parameters.  Note  that  although  the  percent  change  in  damping 
constant  can  te  quite  large,  the  resonant  frequency  changes  no  more  than 
±5%.  This  result  supports  the  proposition  that  first  layer  resonances 
are  associated  with  the  scatterer  itself  not  the  scatterer-ground  plane 
interactions. 

The  real  and  imaginary  part  of  the  natural  modal  current  distribu¬ 
tions  for  resonant  frequencies  y-|i>  y-j ^ >  Y-, 3*  y-j/i  and  y^  are  shown  in 
Figures  3-24  and  3-25.  These  distributions  are,  of  course,  influenced 
by  ground  characteristics  and  scatterer  height-to-length  ratio,  but 
numerical  results  show  that  these  influences  are  relatively  minor.  One 
should  also  observe  that  the  imaginary  part  cf  the  mode  vectors  is  at 
least  an  order  of  magnitude  less  than  the  real  part  indicating  that  the 
mode  function  is  approximately  a  real  function  of  position. 

Coupling  coefficients  associated  with  singularities  y^  y-|2>  and 
y13  are  shown  in  Figures  3-26  through  3-90.  In  each  of  these  figures,  the 
coupling  coefficient,  normalized  such  that  its  maximum  magnitude  is  equal 
to  unity,  is  plotted  as  function  of  the  angle  of  incidence  of  the  electro¬ 
magnetic  excitation,  ie  (the  angle  0  shown  in  Figure  2-1).  In  addition, 
the  curves  arc  presented  with  either  scatterer  height-to-length  ratio,  or 
ground  conductivity,  or  ground  relative  permittivity  as  a  parameter. 
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Figures  3-26  to  3-31  present  coupling  coefficients  for  fundamental 
self-resonant  singularity  y-j -j  at  scatterer  height-to-length  ratios  of 
0.25,  0.50,  0.75,  1.0,  1.25,  and  1.50  respectively.  In  each  of  these 
figures,  the  relative  permittivity,  eR,  is  held  at  unity,  while  the 
curves  in  a  given  figure  correspond  to  ground  conductivities  of  c=120.0, 
a=0.12,  and  o=0. 00012.  When  o=0. 00012  and  eR=1.0,  the  coupling  coef¬ 
ficient  remains  the  same  regardless  of  the  scatterer  height-to-length 
ratio;  an  expected  result  since  the  ground  plane  has  vanished  and  "free 
space"  conditions  prevail.  For  the  "free  space"  case,  the  coupling 
coefficient  has  its  maximum  value  at  broadside  incidence,  ie  (e=90°), 
a  result  which  agrees  with  previous  investigation  [1].  One  observes 
that  the  coupling  coefficient  does  not  differ  greatly  from  the  "free 
space"  case,  Figures  3-26  and  3-27,  even  when  the  ground  conductivity 
is  high,  a=!20.0,  as  long  as  the  scatterer  is  within  approximately  1/4 
wavelength  of  the  ground  plane.  In  Figures  3-28  through  3-31,  the 
coupling  coefficients  for  a=0.12  and  120.0  begin  to  differ  significant' y 
from  the  free  space  case.  With  a  height-to-length  ratio  o*  0.75  and 
a=120.0,  see  Figure  3-28,  maximum  coupling  occurs  at  75°,  and  for  a 
height-to-length  ratio  of  1.00  and  o=120.0,  see  Figure  3-29,  maximum 
coupling  occurs  at  45°.  Now  when  o=120.G  and  the  scatterer  height-to- 
length  ratio  is  1.25  and  1.50  in  Figures  3-30  and  3-31  respectively, 
the  coupling  coefficient  again  maximizes  at  6=90°. 

One  might  logically  ask,  for  a  fixed  scatterer  geometry  and  ground 
parameters,  what  angle  of  incidence  will  produce  maximum  coupling? 
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First  consider  the  simple  case  of  a  time  harmonic  plane  wave  obliquely 
incident  on  a  perfectly  conducting  half-space.  If  the  normal  to  the 
half-space  is  in  the  x-direction,  and  if  the  incident  electric  vector 
is  parallel  to  the  plane  of  incidence,  (vertical  polarization),  then 
the  sum  of  the  incident  and  reflected  waves  will  produce  a  standing 
wave  in  the  x-direction  whose  electric  vector  is  parallel  to  the 
perfectly  conducting  ground  plane.  Defining  a  as  the  angle  formed  by 
the  incident  ray  and  the  ground  plane,  the  magnitude  of  this  standing 
wave  will  be  proportional  to  sin[gx  sina]  sina,  where  e  =  a  the 
free-space  wavelength,  and  x  is  the  perpendicular  distance  from  the 
ground  plane.  Now  let  us  immerse  a  "thin-wire"  in  this  standing  wave 
with  its  axis  parallel  to  the  electric  field  at  a  fixed  number  of 
free-space  wavelengths  above  the  ground  plane.  What  incident  angle  a 
will  produce  a  maximum  standing  wave  at  the  position  of  the  thin  wire? 
It  will  occur  at  the  angle  a  that  maximizes  sin[ex  sina]  sina,  which 
occurs  when  either  cosa=0  or  tan[Bx  sina]  +  ex  sina  =  0. 

Returning  to  the  scatterer-ground  plane  problem  with  the  above 
discussion  in  mind,  will  afford  some  interesting  results.  Table  3-1 
predicts  the  angle  e  that  results  in  maximum  coupling  for  a  given 
scatterer  height  when  the  ground  plane  is  nearly  perfectly  conducting, 
(a=l 20. 0) .  For  Figures  3-26  and  3-27  where  the  scatterer  is  0.11  and 
0.22  wavelengths  above  the  ground  plane,  the  table  predicts  maximum 
coupling  to  occur  at  broadside  incidence  From  Figure  3-28,  maximum 
coupling  occurs  at  e=75°;  for  this  case  the  table  predicts  a  value 
between  73.0°  and  73.5°.  The  difference  in  the  predicted  value  and 


37 


the  value  in  the  figure  is  probably  due  to  numerical  evaluation  of  the 
coupling  coefficient  at  5°  intervals.  Using  the  standing  wave  analogy, 
see  Table  C-l,  one  would  expect  the  coupling  coefficient  to  maximize 
for  42. 0o<u<42. 5°  in  Figure  3-29.  In  the  figure,  the  maximum  occurs 
at  4F°.  Figure  3-30  indicates  peak  coupling  at  0=90°,  a  value  which 
is  predicted  by  the  table.  Although  the  table  also  predicts  a  maximum 
at  approximately  33°,  as  can  be  seen  from  Fig  re  3-30,  there  is  a 
relative  extrema  in  the  magnitude  of  the  coupling  coefficient  at  this 
angle.  Generally,  the  discussion  applicable  to  Figure  3-30  applies  as 
wel 1  to  Figure  3-31 . 

Figures  3-32  and  3-33  display  the  angular  variation  in  the  real 
and  imagine ry  part  of  the  coupling  coefficient  for  the  fundamental 
resonance  with  the  scatterer  height-to-length  ratio  as  a  parameter. 

For  these  figures,  the  conductivity  is  large,  o= 120.0,  and  the  relati\e 
permittivity  is  one,  eR=1.0.  Note  that  although  the  information 
contained  in  these  figures  is  redundant,  (see  Figures  3-26  through  3-11), 
there  inclusion  provides  ready  visualization  of  the  variation  in  coupling 
coefficient  with  scatterer  height-to-length  ratio. 

Figures  3-34  through  3-42  show  variation  in  the  coupling  coefficient 
for  y.| i  with  the  relative  permittivity  as  a  parameter. 

In  Figures  3-34,  3-35  and  3-36,  the  scatterer  height-to-length  ratio 
is  held  at  0.25,  while  the  conductivity  j  is  respectively  120.0,  0.12. 
and  0.00012.  In  each  of  these  figures,  the  coupling  coefficient,  (reel 
and  imaginary  part),  is  plotted  for  a  ground  relative  permittivity,  e^. 
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of  1.0,  15.0,  and  100.0.  As  can  be  seen  from  Figure  3-34,  when  the 
ground  conductivity  is  large,  o=120.0,  there  is  no  variation  in  the 
coupling  coefficient  as  the  ground  permittivity  is  varied.  This 
result  can  be  explained  as  follows.  As  previously  discussed,  the 
singularity  location  relative  to  its  free  space  position  is  a  function 
of  the  discontinuity  at  the  free  space-lossy  ground  plane  interface. 
Although  the  ground  relative  permittivity  for  conditions  depicted  in 
Figure  3-34  varies  greatly,  the  conductivity  remains  large,  producing 
a  large  discontinuity  in  the  ground  plane.  Thus,  the  position  of 
in  the  complex  plane  will  be  constant.  Since  the  coupling  coefficient 
is  strongly  dependant  on  singularity  location,  (see  Equation  2.71),  it 
too  will  be  constant.  In  Figure  3-35,  the  conductivity  o  is  0.12. 
Since  this  is  still  a  relatively  large  value  of  conductivity  only  a 
slight  varation  in  the  coupling  coefficient  is  produced  by  varying 
the  relative  permittivity.  In  Figure  3-36,  the  greatest  variation  in 
the  coupling  coefficient  with  relative  permittivity  is  observed.  The 
conductivity  for  this  case  is  essentially  zero  and,  therefore,  the 
discontinuity  in  the  lossy  ground-free  space  interface  is  controlled 
by  the  relative  permittivity. 

Figures  3-37  through  3-39  and  Figures  3-40  through  3-42  contain 
plots  of  tne  variation  in  coupling  coefficient  with  relative  permit¬ 
tivity  for  scatterer  height- to-length  ratios  of  0.75  and  1.25  respec¬ 
tively.  The  interpretation  of  these  figures  is  similar  to  that  given 
for  Figures  3-34  through  3-36. 
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The  coupling  coefficients  associated  with  the  second  self-resonant 
singularity,  y.^,  are  shown  in  Figures  3-43  through  3-66. 

The  coupling  coefficients  corresponding  to  the  second  sel f-resona  it 
singularity  locations  "A",  "D",  "G"  of  Figure  3-8  are  shown  in  Figures 
3-43  to  3-47.  Table  3-2  predicts  possible  angles  of  maximum  coupling 
for  these  figures  using  the  "standing  wave  analogy".  From  the  first 
row  in  the  table,  it  is  seen  that  there  is  no  angle  e  between  0  and  90 
degrees  that  satisfies  the  equation  TAN[ehSINe]  +  ehSlNo=0.  The  reason 
for  this  is  that  0  <_  ehSINe  £  1.38  for  0°  <  e  <90°  and  thus  TAN[shSINe] 
will  be  positive.  However,  as  always,  a  relative  extrema  is  predicted 
at  e=90°,  an  expected  result  considering  the  fact  that  all  the  coupling 
coefficients  are  even  functions  about  6=90°.  That  is,  one  would  expect 
the  same  amount  of  coupling  at  say  e=95°  as  at  e=85°.  From  Figures 
3-44  and  3-45  maximum  coupling  occurs  at  50  and  35  degrees  respectively, 
values  that  are  predicted  reasonably  well  by  Table  3-2.  Using  the 
"standing  wave  analogy"  one  would  expect  peak  coupling  to  occur  for  an 
incident  angle  of  approximately  22  or  67  degrees  for  the  coupling 
coefficient  in  Figure  3-46.  From  the  figure  peak  coupling  occurs  at 
65°,  with  a  relative  extrema  in  the  magnitude  of  the  coupling  coefficiant 
occuring  at  approximately  22  degrees.  Results  from  the  last  row  of  tha 
table  agree  reasonably  well  with  the  corresponding  figure. 

Variations  in  the  real  and  imaginary  part  of  the  coupling  coefficient 
for  the  second  self-resonant  singularity,  as  a  function  of  spacing, 
h/z,  are  shown  in  Figures  3-48  through  3-51.  In  each  of  these  figures, 
the  relative  permittivity  is  one.  The  conductivity  in  Figures  3-48  and 
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3-49  is  large,  cr=120.0,  while  for  Figures  3-50  and  3-51  the  conductivity 
is  0.12. 

Figures  3-52  through  3-66  are  plots  of  the  coupling  coefficients 
for  singularity  with  the  relative  permittivity  of  the  ground  plane 
as  a  parameter.  Coupling  coefficients  in  these  figures  are  for 
singularity  positions  "A",  "D",  and  "G"  in  Figure  3-8,  and  positions 
"A",  "D",  and  "F"  in  Figure  3-9.  As  before,  when  the  ground  conductivity 
is  large,  the  coupling  coefficient  is  insensitive  to  changes  in  the 
relative  permittivity  (see  Figures  3-52,  3-55,  3-58,  3-61,  and  3-64). 

The  largest  variation  in  the  coupling  coefficient  with  relative  permit¬ 
tivity  is  noticed  in  Figure  3-54,  3-57,  3-60,  3-63,  and  3-66,  a  result 
indicative  of  the  large  difference  in  position  of  points  "G"  and  "F"  in 
Figures  3-8  and  3-9. 

In  Figures  3-67  through  3-90  are  presented  the  coupling  coefficient 
for  the  third  self-resonant  singularity,  y^. 

Figures  3-67  through  3-71  present  the  angular  variation  in  the 
coupling  coefficient  for  singularity  positions  "A",  "D",  and  "G"  in 
Figure  3-10.  Table  3-3  predicts  possible  angles  of  maximum  coupling 
for  these  figures  using  the  "standing  wave  analogy".  From  Figures  3-6? 
and  3-68  maximum  coupling  for  perfect  grcund  (a=120.0)  occurs  at  45 
degrees  and  as  always,  a  relative  extreme  in  the  magnitude  of  the 
coupling  coefficient  occurs  at  90  degrees.  Of  course,  the  extrema  at 
90°  is  predicted  by  Table  3-3,  the  maximum  at  45°,  however,  is  not  since 
BhSINe  <  tt/2  when  h/i  <  0.2.  It  is  also  interesting  to  note  that  unli'<e 
the  corresponding  figures  for  y^  and  yp  in  Figures  3-67  and  3-68  the 
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coupling  coefficient  for  a=0.12  (real  and  imaginary  part)  corresponds 
more  closely  to  the  "free  space"  coefficient  than  to  the  "perfect 
ground"  couoling  coefficient.  This  result  can  be  explained  by  noting 
that  for  higher  order  poles  the  ground  conductivity  must  be  larger 
to  produce  a  perfectly  conducting  earth  than  for  lower  order  poles. 

Next,  consider  Figure  3-69,  for  the  high  conductivity  case  maximum 
coupling  occurs  at  e=40  degrees;  using  the  "standing  wave  analogy" 
one  would  expect  peak  coupling  to  occur  at  approximately  50  degrees. 

The  mechanism  for  this  error  is  not  known.  For  Figures  3-70  and 
3-71  Table  3-3  again  predicts  relatively  accurately  the  angle  of 
maximum  coupling. 

Figures  3-72  through  3-75  show  variation  in  the  real  and  imaginary 
part  of  the  coupling  coefficient  for  the  third  self-resonant  singularity 
with  the  scatterer  height- to-length  ratio  as  a  parameter.  In  Figures 
3-72  and  3-73  the  conductivity  is  120.0  and  the  relative  permittivity 
is  one.  Ir.  Figures  3-74  and  3-75,  the  conductivity  is  0.12  and  the 
relative  permittivity  is  one. 

Figures  3-76  to  3-90  show  variation  in  the  coupling  coefficient 
for  y-ji  with  the  relative  permittivity  as  a  parameter.  In  each  of  these 
figures,  the  conductivity  and  scatterer  height-to-length  ratio  are 
held  constant.  Explanation  of  the  behavior  of  the  coupling  coefficierts 
in  these  figures  is  similar  to  that  given  for  Figures  3-34  through  3-42. 

For  the  wire  scatterer  in  free  space,  the  coupling  coefficient  for 
the  third  singularity  at  an  incident  angle  of  70  degrees  was  shown  by 
Tesche  [1]  to  be  zero.  Note  from  Figures  3-67  through  3-90  that  a  nul 
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in  the  coupling  coefficient  occurs  at  70°  independant  of  the  ground 
parameter  and  scatterer  height- to-length  ratio. 
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Figure  3-2.  Location  of  the  singularities  of  the  thin- 
wire  scatterer  above  a  perfectly  conducting 
ground  plane. 


Figure  3-6.  Trajectory  of  fundamental  self-resonant 

singularity,  y^  ,  as  a  function  of  spacing 

h/£,  and  conductivity,  a.  The  relative 
permittivity  is  held  at  fifteen. 
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Figure  3-12.  Trajectory  of  fourth  sel f- resonant  singularity, 
Y14,  as  a  function  of  spa<  ing,  h /£,  and 

conductivity,  a.  The  relative  permittivity  is 
held  at  one. 


Figure  3-13.  Trajectory  of  fourth  sel f- resonant  singularity 
as  a  function  of  sparing,  h/£,  and 

conductivity,  o.  The  relative  permittivity  is 
held  at  fifteen. 


Figure  3-25.  Real  and  imaginary  part  of  normalized  natural 
modal  current  distributions  for  resonant 
frequencies  y-|4  and  y1  r, -  h/C  =  0.5,  cR  =  1.0, 

a  =  120.0,  0  =  10.6. 


IV.  Conclusions 


The  Singularity  Expansion  Method  (SEM)  coupled  with  the  reflection 
coefficient  approximation  seems  to  be  an  appropriate  technique  for  ana¬ 
lyzing  a  cylindrical  scatterer  over  a  finitely  conducting  ground.  Pole 
locations  of  the  cylinder  over  perfect  ground  calculated  by  investigator 
were  found  to  be  identical  to  those  calculated  in  this  analysis  when  the 
ground  conductivity  is  large.  Also,  when  the  ground  conductivity  is 
reduced  with  the  relative  permittivity  held  at  unity,  calculated  pole 
locations  agree  with  those  of  the  free  space  problem.  These  facts  add 
credence  to  the  reflection  coefficient  technique.  In  general,  it  is 
found  that  displacement  of  the  singularities  from  their  free  space  posi¬ 
tion  is  a  function  of  discontinuity  in  the  ground  plane.  The  mechanism 
for  discontinuity  is  seen  to  be  somewhat  immaterial.  Mode  vectors 
corresponding  to  first  layer  singularities  show  little  or  no  dependence 
on  the  parameters  of  the  problem. 

The  greatest  value  of  this  work  lies  m  the  fact  that  a  relatively 
complicated  problem  has  been  solved  within  acceptable  engineering  accu¬ 
racy  over  a  wide  range  of  parameters. 
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Figure  3-26.  Coupling  coefficient  for  fundamental  self-resonant  singularity, 
Yl i ,  as  a  function  of  conductivity,  j,  and  incident  angle  e. 

The  scatterer  height- to-length  ratio,  h/t,  is  0.25,  and  the 
relative  permittivity,  e„,  is  unity 
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Table  3-1.  Predicted  angle  of  incidence  for  maximum  coupling  using  the 
"standing  wave  analogy. "  The  table  applies  to  Figure  3-26 
through  3-31  for  large  couductivity,  (o  =  120.0),  and  unity 
permittivity,  (eR  =1.0) 
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Imaginary  part  of  coupling  coefficient  for  fundamental  self¬ 
resonant  singularity,  v-| ,  ,  as  a  function  of  scatterer  height-to- 

length  -"ati  o ,  h/ ■  ,  and  incident  angle  The  ground  conductivit 
a,  is  120.0,  and  the  relative  permittivity,  cn,  is  unitv 
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Figure  3-34.  Coupling  coefficient  for  fundamental  self-resonant  singularity, 

y.| i ,  as  a  function  of  the  ground  plane  relative  permittivity,  Ep, 

and  incident  angle  The  scatterer  l,eight-to-length  ratio,  h/£, 
is  0.25,  and  the  conductivity,  j,  is  120.0 
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Figure  3-35.  Coupling  coefficient  for  fundamental  self-resonant  singularity, 

Y-j ^ ,  as  a  function  of  the  ground  pla  e  relative  permittivity,  e^, 

and  incident  angle  o.  The  scatterer  height-to-1 ength  ratio,  h/s., 
is  0.25,  and  the  conductivity,  c,  is  0.12 
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Figure  3-36.  Coupling  coefficient  for  fundamental  self-resonant  singularity, 

Y-j  i ,  as  a  function  of  the  ground  pi  a :  e  relative  permittivity,  e^, 

and  incident  angle  e.  The  scatterer  height-to-length  ratio,  h/t, 
is  0.25,  and  the  conductivity,  c,  is  0.00012 
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Figure  3-38.  Coupling  coefficient  for  fundamental  sel f-resonant  singularity, 
y-i i ,  as  a  function  of  the  ground  plane  relative  permittivity,  c 

and  incident  angle  e.  The  scatterer  heicht-to-length  ratio,  h/ 
is  0.75,  and  the  conductivity,  a,  is  0.12 
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Figure  3.39.  Coupling  coefficient  for  fundamental  sel f-resonant  singular 
y. p  as  a  function  of  the  ground  plane  relative  permit tivit 

and  incident  angle  s.  The  scatterer  height-to-length  ratio 
is  0.75,  and  the  conductivity,  or ,  is  O.CGOlti 


COUPLING  COEFFICIENT 

.00  -0.50  0.00  0.50 


Figure  3-40.  Coupling  coefficient  for  fundamental  self-resonant  singularity, 
y-j i »  as  a  function  of  the  ground  plane  relative  permittivity,  cR 

and  incident  angle  a.  The  scatterer  height-to-length  ratio,  h/i 
is  1.25,  and  the  conductivity,  a,  is  120.0 
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Figure  3-41.  Coupling  coefficient  for  fundamental  sel f-resonanc  sinc.lar 


y 1 1 .  os  a  function  of  the  ground  plane  relative  pern.ittivit 


and  incident  angle  .  Tho  scatterer  height-to-1  ength  ratio 
is  1.25,  and  the  conductivity.  .  is  0.12 
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Figure  3-42.  Coupling  coefficient  for  fundamental  self-resonant  singularity, 
y  1 1 1  as  i  function  of  the  ground  pic-'  j  relative  permittivity,  cR 

and  incident  angle  e.  The  scatterer  ',eicht-to-ler.gth  ratio,  h/?. 
is  1.25,  and  the  conductivity,  c,  is  J. 00012 
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Figure  3-43.  Coupling  coefficient  for  second  self-resonant  singularity,  v]?) 

as  a  function  of  conductivity,  a,  and  incident  angle  •?.  The 
scatterer  height-to-length  ratio,  h/  ,  is  0.20,  and  the  relative 
permittivity,  cR,  is  unity 
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Figure  3-44.  Coupling  coefficient  for  second  sel f- resonant  singularity,  Yp, 

as  a  function  of  conductivity,  o  and  incident  angle  r.  The 
scatterer  heigiit-to-1 ength  ratio,  h/2,  is  0.40,  and  the  relative 
permittivity,  eD,  is  unity 
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Figure  3-47.  Coupling  coefficient  for  second  self- resonant  singularity,  y^, 

as  a  function  of  conductivity,  a,  anc  incident  angle  e.  The 
scatterer  height-to-1 ength  ratio,  h/>  ,  is  1.0,  and  the  relative 
permitt  vity,  cD,  is  unity 
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Table  3-2.  Predicted  angle  of  incidence  for  maximum  coupling  using  the 
"standing  wave  analogy."  The  table  applies  to  Figure  3-43 
through  3-47  for  large  conductivity,  (a  =  120.0),  and  unity 
permittivity,  (e^  =  1.0) 
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Figure  3-43.  Real  part  of  coupling  coefficient  for  second  self-resonant  singu¬ 
larity,  y-|2»  as  a  'Function  of  scatterer  height-to-length  ratio,  h/ ^ 

and  incident  angle  9.  The  ground  cor ductivi ty,  o,  is  120. 0,  and 
the  relative  permittivity,  en,  is  unity 


'0.00  15.03  30.00  45.00  60.00  75.00  90. 

ANGLE  DF.  I  N C  I  D LL N C E 


Figure  3-50.  Real  par;  of  coupling  coefficient  for  second  sel f-resorant  singu¬ 
larity,  ,-| 2 »  as  a  function  of  scatter ?r  height-to-length  ratio,  h./ 

and  incident  angle  9.  The  grouna  conductivity,  a,  is  0.12,  and 
the  relative  permittivity,  eD,  is  unity 
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Figure  3. 51 


Imaginary  part  of  coupling  coeff icier '  for  second  self-resonant 
singular-  ty,  y^,  as  a  function  of  sc-.tterer  height-to- 1  ength  rati  3, 

h/z,  and  incident  angle  e.  The  grourd  conductivity,  c,  is  0.12, 
and  the  relative  permittivity,  e  is  unity 
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Figure  3-52.  Coupling  coefficient  for  second  sel f- resonant  singularity,  -Tj 2 » 

as  a  furction  of  the  ground  plane  re'ative  permittivity,  cD,  and 

K 

incident  angle  9.  The  scatterer  hei<  ht-to-1 ength  ratio,  h/a,  is 
0.20,  erd  the  conductivity,  c,  is  120.0 
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Figure  3-53.  Couplirg  coefficient  for  second  sel f~ 'esonent  singularity,  y^- 

as  a  function  of  the  ground  plane  relative  permittivity,  and 

incident  angle  0.  The  scatterer  heig  it-to-l ength  ratio,  h/t,  is 
0.20,  and  the  conductivity,  a.  is  0.12 
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Figure  3-55.  Coupling  coefficient  for  second  sel f- resonant  singularity, 

as  a  fun:tion  of  the  ground  plane  relative  permittivity,  e^,  and 

incidert  angle  o.  The  scatterer  height- to-length  ratio,  h/t,  is 
0.40,  and  the  conductivity,  a,  is  12C.0 
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Figure  3-56.  Coupling  coefficient  for  second  sel f- resonant  singularity,  y-^, 

as  a  function  of  the  ground  plane  relative  permittivity,  en,  and 

K 

incident  angle  9.  The  scatterer  heic ht-to-1 ength  ratio,  h/t,  is 
0.40,  and  the  conductivity,  c,  is  0.12 
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Figure  3-57.  Coupling  coefficient  for  second  self-resonant  singularity, 

as  a  function  of  the  ground  plane  relative  permittivity,  cR 

incident  angle  6.  The  scatterer  height-to-length  rat':n,  h/ 
0.40,  and  the  conductivity,  c,  is  O.C0012 
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Figure  3-50.  Couplirg  coefficient  for  second  sel f-resonant  singularity, 

as  a  function  of  the  ground  plane  relative  permittivity,  e 

incident  angle  3.  The  scatterer  height- to-length  ratio,  h/ 
0.60,  and  the  conductivity,  c,  is  120.0 


COUPLING  COEFFICIENT 
.00  -0.50  0.00  0.50 


COUPLING  COEFFICIENT 

.00  -0.50  0.00  0.50 


«b=1.0,RERI.-3  ,  IMfiG.-A' 
SR=lS.O,REni-+  .IKRG.-X 


60.00 


0.00  15. CO  30.00  45.00  60.  i 

RNGLE  OF  INCIDENCE 


75.00 


90.00 


Figure  3-60.  Coupling  coefficient  for  second  self-resonant  singularity, 

as  a  function  of  the  ground  plane  relative  permittivity,  and 

incident  angle  e.  The  scatterer  heicht-to-length  ratio,  h/£,  is 
0.60,  c.nd  the  conductivity,  o,  is  0.00012 
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Figure  3-62.  Coupling  coefficient  for  second  self-resonant  singularity 

as  a  function  of  the  ground  plane  relative  permittivity, 

incident  angle  e.  The  scatterer  heignt-to-length  ratio, 
0.80,  and  the  conductivity,  o,  is  0.12 
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gure  3-63.  Coupling  coefficient  for  second  sel f-resonant  singularity,  y 

as  a  function  of  the  ground  plane  relative  permittivity,  c^, 

incident  angle  o.  The  scatterer  hei ;  ht-to-lenqth  ratio,  h/ 
0.30,  end  the  conductivity,  o,  is  0.00012 


COUPLING  COEFFICIENT 

.00  -0.50  0.00  0.50 


COUPLING  COEFFICIENT 

.CO  -0.50  0.00  0.5C 


«r*1.0,RERL-O  . ItfHG. -A 
€a=  1 5 . 0 ,  REAL--f  , I^AG.-x 


0.00  15.00  30.00  U5.00  60.00 

ANGLE  OF  INCIDENCE 


75.00 


Figure  3-65.  Couplirg  coefficient  for  second  self -resonant  singularity, 
as  a  function  of  the  ground  plane  relative  permittivity,  eD 

K 

incident  angle  9.  The  scatterer  hei:ht-to-length  ratio,  h/ 
1.0,  and  the  conductivity,  a,  is  0.1 7 
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Figure  3-56.  Coupling  coefficient  for  second  self-resonant  singularity, 

as  a  function  of  the  ground  olane  relative  oermittivity,  en,  and 

J  R 

incident  angle  e.  The  scatterer  height-to-1 ength  ratio,  h/t,  is 
1.0,  and  the  conductivity,  c,  is  0.0C012 
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gure  3-67.  Coupling  coefficient  for  third  sel1"-  -esor.ant  singularity,  as 

a  function  of  conductivity,  n  and  in  ident  angle  .  The  scatter? 
height-fo-1 ength  ratio,  h/t,  is  0.10  and  tne  relative  permittivi 
eD ,  is  unity 
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Figure  3-71.  Couplinc  coefficient  for  third  sel f-- esonant  singularity,  y-, 

a  function  of  conductivity,  o',  and  i  cident  :ncle  e.  The  sc 
height-to-length  ratio,  h/t,  is  0.50,  and  the  relative  p c-rtrii 
e0,  is  unity 
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Table  3-3.  Predicted  angle  of  incidence  maxima  :oupling  t'-e 

("standing  wave  analog;  .  The  table  applies  to  F-  .res  2-z~ 

through  3-71  for  larg:  ccrductiv ■ ty,  (o  -  12U.0),  and  „n':_ 
permittivity,  (tp  =  i.O) 
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Figure  3-72.  Real  part  of  coupling  coefficient  f  third  sel F-i  it  singul 
y  1 3  •»  as  •  ■  - .  -  - '  io,  h/i,  nd 

incident  angle  '.  The  grout'.'  conductivity,  c,  is  120.0,  and  the 
relative  penni :ti.i ty ,  c,,,  is  unity 
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Figure  3-74.  Real  pe^l.  of  coupling  coefficient  fc-'  third  jelf-resor.ant  singularity, 
v13,  as  c  function  of  scatterer  hMgh  -tc-length  ratio,  h /z,  and 

incident  angle  e.  The  ground  cor  uct  vity,  a,  is  0.12,  and  the 
relative  permittivity,  ed,  is  unity 
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76.  Coupling  coefficient  for  third  self-resonant  since! 

a  function  of  the  ground  plane  relative  permittivity 

angle  e.  The  scatterer  hsight-to-1  igth  ratio,  h/  , 
the  concucti vity,  c,  is  120.0 
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Figure  3-71 
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Coupling  coefficient  for  third  sel f- resonant  singularity,  y.,,  a; 

a  functor  o.  the  ground  plane  relat've  permittivity,  en,  and  inc  dent 

angle  e.  The  scatte-er  height-to-length  ratio,  h/t,  is  0.10,  end 
the  _ oncuctivity,  ,  is  0.12 


78.  Coupling  coefficient  for  third  self-resonant  singularity, 
as  a  function  ot  the  ground  plane  relative  permittivity,  r 

H 

incident  angle  .  The  scatterer  heicht-to-length  ratio,  h/ 
0.10,  and  the  conductivity,  is  0.C0012 
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Figure  3-79.  Coupling  coefficient  for  third  sel f- resonant  singularity, 

as  a  function  of  the  ground  plane  relative  permittivity,  e 

incident  angle  o.  The  scatterer  heic ht-to-length  ratio.  It 
0.20,  ard  the  conductiv i ty ,  o,  is  12C.0 
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Figure  3-80.  Couplinj  coefficient  for  third  self -resoni  t  singularity  v  , 

as  a  function  of  the  ground  plane  relative  p:< v.ittwity,  eR,  and 

inciden  angle  e.  The  scatter er  hei  -ht-to-1 ength  ratio,  h/e,  is 
0.20,  and  the  conductivity,  o,  is  0.  2 
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Figure  3-81.  Couplin:  coefficient  fot  third  self-resonant  singularity,  y]  v 

as  a  function  of  the  ground  plane  re  ative  permittivity,  eR,  and 

incident  angle  9.  T.,e  scatterer  heicht-to-length  ratio,  h/£,  is 
0.20,  aid  the  conductivity,  a,  is  0.1'C012 


Figure  3-82.  Couplin:  coefficient  for  third  sel f- resonant  singularity,  v-p, 

as  a  function  of  the  ground  plane  relative  permittivity,  ,  and 

incident  angle  e.  The  scatterer  he i c ht-to-length  ratio,  h/i,  is 
0.30,  a'd  the  conductivity,  a,  is  120.0 
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Figure  3-83. 


Coupling  coefficient  for  third  sel f- 'esonant  singularity,  y 

as  a  function  of  the  ground  plane  relative  permittivity,  ep,  and 

incident  angle  6.  The  scatterer  height-to-ler.gth  ratio,  h/a,  is 
0.30,  ani  the  conductivity,  a,  is  0.12 
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Figure  3-84.  Couplin:  coefficient  for  third  sel f- resonant  singularity, 
as  a  function  of  the  ground  plane  relative  permittivity,  cD 
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incident  angle  r .  The  scatterer  heic ht-to-1 enqth  ratio,  h/ 
0.30,  a-d  the  conductivity,  c,  is  0.(0012 
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figure  3-89.  Coupling  coefficient  for  third  seif-resonant 


as  a  unction  of  the  ground  plane  'dative  per 
h/£,  is  0.50,  and  the  conductivity  o,  is  0.18 
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Figure  3-30.  Coup!  ng  coefficient  for  third  se  f-resonant  singularity 
as  a  unction  of  the  ground  plane  'elative  pernittivi ty, 
h/t,  s  0.50,  and  the  conductivity  a ,  is  0.00012 


IV.  Cor  el  us i or  3 


The  Singularity  Expansion  Method  (SEM)  coupled  with  the  reflection 
coefficient  approximation  seers  to  be  an  appropri  te  -•  chni  le  for  ana¬ 
lyzing  a  cylindrical  scatterer  over  a  firitely  conducing  ground.  Pole 
locations  0*  the  cylinder  yrer  perfect  ground  calculated  in  Ligators 
were  found  :o  be  identical  to  those  c;lculated  in  this  analysis  when  the 
ground  conductivity  is  large.  Also,  wher  the  ground  conductivity  is 
reduced  witu  the  relative  permittivity  held  at  unity,  calculated  pole 
leer  ti  s  agree  with  these  of  the  free  space  problem.  These  facts  add 
credence  to  the  reflection  coefficient  technique.  In  general,  it  is 
found  that  displacement  of  the  singularities  from  their  free  space  pos¬ 
tion  is  a  f  notion  of  discontinuity  in  tie  ground  plane.  The  mechanism 
for  discontinuity  is  seen  to  be  somewhat  immaterial.  Mode  vectors 
corresponding  to  first  layer  singularities  show  little  or  no  dependence 
on  the  parameters  of  the  problem. 

Coupling  coefficients  corresponding  a  the  first  three  fundaments 
resonances  lave  been  presented  as  a  function  of  the  parameters  of  the 
thin  wire-lossy  ground  problem.  Tor  the  case  when  the  ground  plane  is 
nearly  perfectly  conducting  the  "standirg  wave  analogy"  as  discussed  in 
the  numeric  1  results  section  predicts,  relatively  accurately,  the  ang  e 
of  incidenc:  which  results  in  ma: imum  coupling  to  the  scatterer.  Although 
the  "standing  wave  analogy"  is  obviously  not  completely  rigorous,  the 
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point  of  view  taken  here  is  that  the  physical  insight  provided  this 
technique  and  its  apparent  accuracy  in  predicting  maximt,  coupling  are 
sufficient  *o  compensate  for  the  1  ck  of  1  ri  -  . 

The  greatest  value  of  this  work  lies  in  the  fact  that  a  relativel, 
cc  plicated  problem  has  been  solved  withi  acce  table  €  >ring 
racy  over  a  wide  range  of  parameters. 
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